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Abstract 

The advent of the Auger Engineering Radio Array (AERA) necessitates the development of a powerful framework for the 
analysis of radio measurements of cosmic ray air showers. As AERA performs "radio-hybrid" measurements of air shower radio 
emission in coincidence with the surface particle detectors and fluorescence telescopes of the Pierre Auger Observatory, the radio 
analysis functionality had to be incorporated in the existing hybrid analysis solutions for fluoresence and surface detector data. 
This goal has been achieved in a natural way by extending the existing Auger Offline software framework with radio functionaUty. 
In this article, we lay out the design, highlights and features of the radio extension implemented in the Auger Offline framework. 
Its functionality has achieved a high degree of sophistication and offers advanced features such as vectorial reconstruction of the 
electric field, advanced signal processing algorithms, a transparent and efficient handling of FFTs, a very detailed simulation of 
detector effects, and the read-in of multiple data formats including data from various radio simulation codes. The source code of 
this radio functionality can be made available to interested parties on request. 

Keywords: cosmic rays, radio detection, analysis software, detector simulation 



1. Introduction 

Forty years after the initial discovery of radio emission from 
extensive air showers Ol, the CODALEMA H and LOPES fl 
experiments have re-ignited very active research activities in the 
field of radio detection of cosmic ray air showers. Nowadays, 
the field is in a phase of transition from first-generation experi- 
ments covering an area of less than 0. 1 km^ to large-scale arrays 
of tens of km^. In particular, the Auger Engineering Radio Ar- 
ray (AERA) 0] will complement the southern site of the Pierre 
Auger Observatory f3\ with 161 autonomous radio detector sta- 
tions covering an area of » 20 km^. 

One particular merit of the Pierre Auger Observatory is its 
hybrid mode of observation, which uses coincident detection of 



extensive air showers with both optical fluorescence telescopes 
(FD) and surface particle detectors (SD) to gain in-depth infor- 
mation on the measured air showers. Consequently, the analysis 
software has to support complete hybrid processing and inter- 
pretation of the data. This requirement is fulfilled by the Auger 
Offline software framework fa]. To take full advantage of the 
radio data taken in the hybrid environment of the Pierre Auger 
Observatory, it is clear that also radio analysis functionality, 
which has so far been existing in a separate software package 
]0], had to be included in this hybrid analysis framework. 

In this article, we describe how we have therefore built ad- 
vanced radio analysis functionality into the Auger Offline soft- 
ware framework. The general structure of the radio implemen- 
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tation in the Offline framework will be discussed in section |2] 
A number of innovative features have been realized in this con- 
text for the very first time. These and other highlights will be 
discussed in section [3] Finally, in section |4] we demonstrate 
how the advanced radio functionality embedded in the Offline 
framework can be used to carry out a complete detector simula- 
tion and event reconstruction on the basis of a simulated radio 
event. 

2. Embedding radio functionality in the Offline framework 

The Offline framework has a clear structure to allow for easy 
maintenance and ongoing shared development over the whole 
Ufe-time of the Pierre Auger Observatory f^. In particular, 
there is a clear separation between the internal representation 
of the Detector and the Event. The Detector provides access 
to all of the relevant detector information such as the positions 
of detector stations in the field, the hardware associated with 
these stations, etc. The Event data structures in contrast hold 
all of the data applying to a specific event, such as ADC traces, 
but also reconstructed quantities such as the event geometry. 
There is no direct connection between these two entities. In- 
stead, analysis Modules use the defined interfaces of both the 
Detector and Event data structures to carry out their specific 
analysis tasks. No interface exists either between separate anal- 
ysis modules, which can only propagate their results through 
the Event data structure. This ensures that dependencies be- 
tween analysis modules are kept to a minimum and facilitates 
the replacement of individual modules with alternative imple- 
mentations, thereby providing a very high degree of flexibility. 

Clearly, the radio analysis functionality had to be imple- 
mented following the same philosophy. The hierarchical im- 
plementation of the radio parts of both the Detector and Event 
classes in addition to the existing FD- and SD-specific classes 
is depicted in Fig.[T] In analogy to the hierarchy of Stations and 
PMTs in the SD functionality, the implementation of the radio 
data structures has been divided into Stations and Channels. A 
Station represents one location in the field at which the electric 
field of the radio waves is measured. Data stored at Station level 
therefore represents the physical electric field devoid of any de- 
tector influence except for the location (and limited observing 
bandwidth) of the Station. In contrast. Channels represent the 
individual antenna channels at which the "raw" measurement is 
performed by an ADC digitizing voltages. This clear separation 
between Channels and Stations is a very powerful concept and 
is original to the radio implementation in Offline. We will dis- 
cuss its significance, among other highlights, in the following 
section. 

3. Highlights of the radio analysis functionality 

The radio functionality in the Offline framework provides a 
number of unique features facilitating an advanced radio data 
analysis. In this section, we will describe some of these high- 
lights. 



3.1. Clear separation of Channel- and Station-levels 

When analyzing radio data, one is faced with two different 
"levels". The Channel level is defined by the detector chan- 
nels acquiring the raw data. These data consist of time-series 
of samples digitized with a sampling rate adequate for the fre- 
quency window of interest. Each sample denotes a scalar quan- 
tity such as an ADC count recorded by the channel ADC. Low- 
level detector effects such as the correction for the frequency- 
dependent response of cables, filters and amplifiers are treated 
on this level for each Channel individually. Likewise, detector- 
related studies such as the evaluation of trigger efficiencies 
would be typically performed on Channel level. When read- 
ing in measured data files, the raw data (ADC counts) are filled 
into the appropriate Channel data structures. 

In contrast, the Station level is defined by the physical elec- 
tric field present at a given location in the field, stored as a 
time-series of three-dimensional vectors. It is on Station level 
that radio pulses are identified and quantified, before a geom- 
etry reconstruction of the given event is performed. Once the 
event reconstruction has been completed, the data at Station 
level no longer have any dependence on the detector character- 
istics, except for the location and limited observing bandwidth 
of the measurement. A reconstruction of the electric field on the 
Station level is therefore suited best for a comparison of radio 
measurements of different experiments, as well as for the com- 
parison of radio measurements with corresponding simulations. 
Since simulated electric field traces provided by radio emission 
models also represent physical electric fields independent of a 
given detector, they are read in on the Station level. 

Analysis modules in Offline usually work on either Channel 
or Station level, and typically it is very clear which analysis step 
has to be performed on which level. The transition between the 
two levels is performed by applying the characteristics of the 
antennas associated to each of the Channels. This transition 
can be employed in both directions, from Station to Channel 
or vice-versa. The transition from a Station to the associated 
Channels is typically performed to calculate the response of 
the individual detector Channels to an electric field provided 
by simulations. The opposite transition is required when re- 
constructing the three-dimensional electric field vector from the 
data recorded by the (typically) two measurement channels in 
the field. This reconstruction will be further discussed in sec- 
tion [O 

3.2. Read-in from different data sources 

The Event data structures are complemented with reader 
functionality to populate them with data available in one of 
several file formats for both experimentally measured data and 
simulated radio event data. Due to its wealth of supported for- 
mats and the possibility of easy extension with new formats, 
the radio functionality in Offline therefore provides very power- 
ful functionality to compare data and simulations from different 
sources, which again is an original feature usually not found in 
the analysis software suites developed in the contexts of other 
experiments. At the time of writing, the following data formats 
are supported. For experimental data: 
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• measurement data from two different prototype setups sit- 
uated at the Balloon Launching Station of the Pierre Auger 
Observatory flU 

• measurement data from AERA 

For simulation data, the following formats are currently read- 
able: 

• simulation data from MGMRlfloll 



simulation data from REAS2 and REAS3 ||11U12|1 



• simulation data from ReAIRES 113 1 



3.3. Modular approach 

The strict interface design of the Detector, the Event and the 
analysis modules allows for a very modular implementation of 
radio analysis functionality. As the analysis modules are the 
part of the code typically the most exposed to the end-user, their 
interface has been kept relatively simple. End-users developing 
analysis functionality for Offline therefore only need relatively 
basic proficiency in C++. 

An analysis application within Offline is defined through a 
"module sequence" in XML syntax, an example of which is 
listed in section |4] In such a module sequence, analysis mod- 
ules are chained in a meaningful sequence to perform a specific 
analysis task. The individual modules do not communicate di- 
rectly with each other, but only share data through the Event 
data structures. Consequently, modules can easily be removed, 
replaced or rearranged within a module sequence. This does not 
require recompilation of the source code. Additionally, each 
module can be configured individually through XML files. 

3.4. Transparent FFT handling 

Radio analyses typically apply algorithms both on time- and 
frequency-domain data. As a consequence, they heavily rely 
on fast Fourier transforms (FFTs). The Offline framework has 
thus been extended with FFT functionality based on the FFTW 
library 11411 . A special feature of this implementation is that 
FFTs are handled completely transparently in the background. 
The user does not need to invoke FFTs manually. 

This is realized by the use of FFTDataContainers as illus- 
trated in Fig.|2] These containers encapsulate both the time- and 
frequency-domain representations of radio data on the Channel 
and Station levels. The user can access both the time-domain 
and frequency-domain data at any time. The FFTDataCon- 
tainer keeps track of which representation has been changed 
last and whether an FFT has to be performed or not before the 
data requested by the user are returned. All data are passed 
by reference and changed in place, so that even traces with an 
extreme length can be handled efficiently. 

As a consequence of this design, the user can simply chain 
analysis modules working in any of the two domains without 
worrying which domain has last been worked on. (There is 
a performance benefit when grouping modules working in the 
same domain together, but it is not very significant.) 



3.5. Advanced analysis modules 

A number of analysis modules performing recurring steps in 
advanced radio analysis pipelines are available by default. They 
can easily be included or excluded from module sequences as 
needed: 

• modules applying bandpass filters to the Channel and Sta- 
tion levels 

• a module performing an up-sampling of under-sampled 
data 

• a module resampling data to a diflterent time-base 

• a module suppressing narrow-band radio frequency inter- 
ference through a "median filter" 

• a module performing an enveloping of time traces via a 
Hilbert transform 

• a module determining timing differences between differ- 
ent antenna stations from the reference phases of a beacon 
transmitter 

• modules applying a windowing function (e.g., Hann win- 
dow) to the Channel and Station levels 

3.6. Detailed simulation of the detector response 

When comparing measured data to simulated radio pulses 
from various models, it is required to perform a detailed simu- 
lation of the effects introduced by the various detector compo- 
nents. This encompasses in particular: 

• the complex response (impulse response defined by the 
frequency-dependent amplitudes and phased of all the 
analogue components (cables, filters, amplifiers) in each 
individual channel 

• the frequency- and direction-dependent complex gain (or 
"effective antenna height") of the antenna connected to 
each individual channel (cf. Fig. [3]) 

• effects introduced by the sampling of the data with a given 
sampling rate 

• saturation effects occurring at the ADCs 

• effects introduced by the layout of the array, including ge- 
ometric effects occurring on large scales due to the curva- 
ture of the Earth 

All of this functionality has been implemented in the Offline 
framework. At the moment, detector description data are pro- 
vided as XML files. Later, a transition to MySQL or SQLite 



full transport matrix representing the transmission in forward direction, 
transmission in backward direction, as well as the reflections on the input and 
output could be implemented for a more detailed description. For the moment, 
however, we assume that impedance matching in the experimental setup is suffi- 
ciently good so that transmission in the forward direction describes the detector 
response with good precision. 
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RecEvent 



Figure 3: The antenna chai'acteristics are defined in a splierical coordinate sys- 
tem with the antenna in its center The effective antenna height fi^ff for emis- 
sion coming from a given arrival direction is decomposed into the components 
along the unit vectors eg and e^. These local (i.e., arrival direction-dependent) 
unit vectors lie in the plane perpendicular to the Poynting vector, which aligns 
with the dashed blue line in this figure for the incoming direction defined by 
and ip. As electromagnetic waves in air have no electric field component along 
the Poynting vector, this representation is complete. The antenna height itself is 
a frequency-dependent, complex quantity, the amplitude of which denotes the 
gain of the antenna, while the phase provides information about signal delays 
and dispersion. 



databases is foreseen and can be performed transparently. The 
complex response of individual Channels is provided via a Re- 
sponseMap detailing the hardware elements comprising each 
individual channel. The overall response of each channel is 
then calculated on-the-fly from the tabulated responses of each 
individual hardware component listed in the ResponseMap. A 
caching mechanism ensures that overall responses are only re- 
calculated when needed. 

3.7. Vectorial E-field reconstruction 

The physical electric field is a three-dimensional, vectorial 
quantity. When comparing results from different experiments 
or experimental results and radio emission models, the electric 
field is the quantity of choice, as in principle it has no depen- 
dence on the detector (except for the location at which it was 
measured and the limited observing bandwidth). Most radio de- 
tectors, however, are only equipped with two channels per po- 
sition in the field, typically measuring the east-west and north- 
south linear polarization components. In other words, they only 
measure a projection of the three-dimensional electric field to 
the horizontal plane. In such a setup, two Channels are avail- 
able at each detector station, one connected with an east-west- 
aligned antenna and one connected with a north-south-aligned 
antenna. The (scalar) response of each individual channel to the 
incoming electric field can be calculated as the dot product be- 
tween the electric field vector and the vectorial representation 
of the effective antenna height. The effective antenna height is 
related to the gain of the antenna and depends on the arrival 
direction and frequency (cf. Fig. |3]l. Consequently, this calcu- 
lation is best done in the frequency domain. 
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Figure 4: Hierarchical structure of the radio-related quantities stored in ADST 
files. The individual classes constitute logical entities present in the Event data 
structures, storing for example information on Monte Carlo generated quanti- 
ties in RGenStation and the reconstructed radio lateral distribution function in 
RdLDF. 



The more difficult problem is the inverse calculation: the re- 
construction of the three-dimensional electric field vector from 
the two-dimensional measurement. This inversion is possible if 
the arrival direction of the electromagnetic wave is known, be- 
cause electromagnetic waves in the atmosphere constitute trans- 
verse waves, the electric field of which lies in a plane perpen- 
dicular to the direction of propagation^ 

The antenna characteristics needed to reconstruct the vecto- 
rial electric field on the Station level depend on the arrival di- 
rection. This arrival direction, however, is not available until 
after the reconstruction on the basis of the Station level three- 
dimensional electric field. Therefore, an iterative approach 
starting with a reasonable initial arrival direction is performed 
in the radio analysis in Offline. The reconstructed arrival direc- 
tion quickly converges to its final value, and the vertical com- 
ponent of the electric field can be reconstructed from the two- 
dimensional measurement (cf. section |4]i. This reconstruction 
scheme is truly original to the Offline radio functionality. 

3.8. Advanced data output 

End-users performing a higher-level analysis usually do not 
need access to all the raw data of each individual event. They 
rather need information about reconstructed quantities derived 
by the low-level analysis pipeline. A data format to store such 
quantities reconstructed by an Offline-based analysis is the Ad- 
vanced Data Summary Tree (ADST). ADST files hold relevant 
data for FD, SD and now also radio-reconstructed quantities. 
The structure of the radio data within the ADST files is de- 
picted in Fig. |4] In addition to accessing the content of ADST 
files directly from end-user analysis programs, a graphical user 
interface exists for browsing the contents of ADST files and vi- 
sualizing the included events. This EventBrowser has also been 
complemented with radio-specific functionality, so that also the 



'*At the moment, we assume that the propagation direction of radio emission 
from extensive air showers can be approximated well with the shower axis. 
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<module> EventFileReaderDG </module> 
<module> RdSt at i onAs soc i at or </module> 

<! -- cf . fig. 6 raw simulated electric fields --> 

<module> RdAnt ennaSt at ionToChanne IConvert er </module> 
<module> RdChannelNoiseGenerator </module> 

<! -- cf . Fig. 7 voltages at antenna foot-points --> 



<module> RdChannelResponselncorporator 
<!-- cf . Fig. 8 voltages at ADC 



</module > 



<module> RdChannelResampler </module> 
<module> RdChanne ITime Ser i e sCl ipper </module> 
<module> RdChannelVoltageToADCConverter </module> 

<! -- signal is now equivalent to raw ADC values --> 



<module> RdChannelADCToVoltageConverter </module> 

<module> RdChannelPedestalRemover </module> 

<module> RdChannelResponselncorporator </module> 

<module> RdChannelRFISuppressor </module> 

<module> RdChannelUpsampler </module> 

<module> RdChannelBandpassFilter </module> 



< ! 



cf . Fig. 9 voltages after signal processing 



<loop numTimes =" unbounded" > 

<niodule> RdAnt ennaChannelToSt at ionConvert er </module> 

<module> RdS t at i onS i gnalRe c onst ruct or </module> 

<module> RdD ir e c t i onConver genceChe cker </module> 

<module> RdPlaneFit </module> 

</loop> 

<module> RdSt at i onWindowSe t t er </module> 
<! -- cf . Fig. 10 field strengths in clipped trace --> 

<module> RdStationTimeSeriesWindower </module> 

<! -- cf . Fig. 11 Hann - windowed field strengths --> 



<module> RecDataWriter 



</module > 



radio part of the event (such as traces, spectra, a map of the ar- , 
ray and event geometry, ...) can be visualized in a user-friendly ^ 
and intuitive fashion. 4 

5 
6 

3.9. Easy extension to hybrid analysis ' 

8 

One of the main motivations for including radio functional- 
ity in the Offline framework was to exploit the hybrid nature of h 
the data acquired within the Pierre Auger Observatory. Han- " 
dling all detector data within the same analysis software will w 
make it easily possible to develop analysis procedures combin- \l 
ing data from SD and radio, FD and radio or SD, FD and radio n 
altogether. This will be possible both for measured data and |* 
for simulated events. Developing the "radio-hybrid" analysis 20 
strategies is yet a challenge for the future. The technical pre- 
requisites for this development have, however, been success- 23 
fully provided with the inclusion of radio functionality in the 
Offline framework. 26 

27 
28 
29 

4. Example for an analysis pipeline ^ 

32 

To illustrate the flexibility and level of sophistication 
achieved with the radio functionality in Offline, we discuss 35 
here a module sequence processing simulated data with a 
reconstruction pipeline that incorporates all relevant detector 38 
effects. Each line in the XML file listed in Fig. |5] denotes 
a module being invoked to perform a specific analysis step. 
Radio modules starting with RdChannel manipulate the event 
data on the Channel level, modules starting with RdStation 
manipulate data on the Station level. The modules starting 
with RdAntenna perform the transitions between the two levels 
(cf. section [TT] ). Modules can easily be removed, replaced or 
rearranged to change the analysis procedure without having to 
recompile the source code. In the following, we will briefly 
discuss the individual steps of the module sequence and show 
how the signal evolves on its way through the analysis pipeline. 



4.1. Read-in and association 

The module sequence starts with the read-in of simulated 
event dat^ using the EventFileReaderDG module. After 
read-in, the simulated data represent an "abstract" simulation 
that is not yet associated to any detector stations. This as- 
sociation is performed by the RdStationAssociatior mod- 
ule, which associates the simulated signal traces with the corre- 
sponding stations in the field, and at the same time pads them 
appropriately to ensure that the signal falls into the correct part 
of the time series trace. Afterwards, the Station data structure 
contains the physical electric field vector as predicted by the 
simulation, without the inclusion of any detector effects. The 
corresponding traces and spectra are depicted in Fig. |6] 



'The example event used here has an energy of 2.1 • lO'* eV, a zenith angle 
of 58.4°, an azimuth angle of 291.0° con'esponding to an anival direction of 
approximately SSW, and has been simulated with a proton as primary particle. 



Figure 5 : Example module sequence for performing a full detector simulation 
and event reconstruction of a simulated radio event. The time series data and 
frequency spectra extracted at certain steps during this analysis pipeline are 
illustrated with the figures referenced in the sequence. Analysis modules can 
be re-ordered or exchanged without any recompilation of the source code. 

4.2. Simulation of the detector response 

The next steps in the module sequence change the data 
such that they become equivalent to data measured experimen- 
tally. The RdAntennaStationToChctnrtelConverter calcu- 
lates the signal voltages that each Channel of a given Station 
would have seen at the foot-points of the corresponding an- 
tennas by folding in the antenna response applicable to each 
individual channel. In the typical case of two antenna Chan- 
nels per Station, this means that the three-dimensional elec- 
tric field vector is projected to a two-dimensional surface. The 
RdChannelNoiseGenerator module then adds broad-band 
radio noise to the event. The resulting, simulated data for the 
east and north channels are shown in Fig.|7] 

The following call of the 

RdChannelResponselncorporator incorporates the (for- 
ward) detector response of the cables, filters and amplifiers 
comprising the corresponding Channel. After this module, the 
signal represents the voltages that would be measured at the 
channel ADCs, depicted in Fig. [8] 

The following steps convert this voltage at the ADCs to 
the signal that the channel ADC would indeed have mea- 
sured. The RdChannelResampler module re-samples the 
Channel time series data to the time-base with which the 
data are sampled in the experiment. (The prerequisite to 
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this module is that high-frequency components which could 
lead to aliasing effects have been suppressed. This is en- 
sured here because the RdChannelResponselncorporator 
includes filters that do just that in the experimental setup.) The 
RdChannelTimeSeriesClipper then clips the Channel traces 
to the number of samples which are taken in the experimental 
setup. Finally, the RdChannelVoltageToADCConverter con- 
verts the voltages of each sample to ADC counts that would 
have been recorded by the channel ADCs, thereby taking into 
account quantization and saturation effects. 

At this point, all of the relevant detector effects have been in- 
corporated in the simulated traces. In other words, the data now 
have the same properties as measured data directly after read- 
in. Consequently, the remainder of the reconstruction pipeline 
is identical to one that would be applied directly to measured 
data. Choosing the reconstruction procedure applied to the sim- 
ulated data identical to the one applied to measured data makes 
sure that even subtle changes introduced by individual analysis 
modules can be investigated on the basis of simulations. 

4.3. Signal cleaning 

After converting the ADC counts back to voltages 
with the RdChannelADCToVoltageConverter and re- 
moving a possible DC offset of the ADC with the 
RdChannelPedestalRemover, the characteristics of the ana- 
logue components of each Channel are folded out from the data 
with a second call of the RdChannelResponselncorporator. 

The following steps are intended to improve the reconstruc- 
tion quality by the use of advanced digital processing tech- 
niques. In a first step, the RdChannelRFISuppressor mod- 
ule suppresses narrow-band signals (e.g., TV carriers present 
in measured data) using a median filter It is followed by the 
RdChannelUpsampler module which correctly reconstructs 
(interpolates) the signal on a finer time-base. (This is possi- 
ble if the complete signal information is present in the digi- 
tized signal, i.e., if the Nyquist criterion for data sampling is 
fulfilled.) In a further step the signal bandwidth is limited by a 
digital bandpass filter using the RdChannelBandpassFilter 
module. After these steps, the signal is ready for the reconstruc- 
tion of the three-dimensional electric field vector and looks like 
the data presented in Fig.|9] 

4.4. Vectorial reconstruction 

The following loop performs an iterative reconstruction of 
the three-dimensional electric field vector and the signal arrival 
direction. The RdAntennaChannelToStationConverter 
performs the reconstruction of the three-dimensional elec- 
tric field described in section 13.71 Afterwards, the 
RdStationSignalReconstructor identifies the times at 
which radio pulses have been detected and quantifies the param- 
eters of these pulses. Next, the RdPlaneFit reconstructs the 
arrival direction of the radio signal with a plane-wave assump- 
tion based on the previously established pulse arrival times. 
Finally, the RdDirectionConvergenceChecker module tests 
whether the iterative procedure has converged or not and breaks 
the loop accordingly. 



4.5. Post-processing 

After breaking the iterative reconstruction loop, the vecto- 
rial, detector-independent electric field has been completely re- 
constructed. For practical purposes, the time series is then 
restricted to a window of 500 ns around the detected pulses 
with the RdStationWindowSetter. This leads to leak- 
age effects in the frequency spectra, visible in Fig. [10] To 
suppress these leakage artifacts, a Hann window is applied 
with the RdStationTimeSeriesWindower. The final recon- 
structed signal is then seen in Fig. [TT] As a last step, the 
RecDataWriter call writes out the reconstructed event data to 
an ADST file for further processing in higher-level analyses. 
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5. Conclusions and Outlook 

We have implemented a complete set of radio analysis func- 
tionality in the Offline software framework of the Pierre Auger 
Observatory. The radio functionality has been included in a 
canonical and seamless way in addition to the existing SD and 
FD functionality. This approach will make the realization of 
"radio-hybrid" analysis strategies in the future straight-forward. 

Already now, however, the radio functionality in Offline has 
reached a high degree of sophistication with highlights such as 
a very fine-grained simulation of detector effects, advanced sig- 
nal processing algorithms, transparent and efficient handling of 
FFTs, read-in of multiple file formats for measured and sim- 
ulated radio data, and in particular the reconstruction of the 
three-dimensional electric field vector from two-dimensional 
measurements. Planned improvements encompass the imple- 
mentation of a curved fit, inclusion of interferometric radio 
analysis functionality, and the handling of a time-variable de- 
tector including a fine-grained treatment of the instrumental cal- 
ibration. 

Parties interested in using the functionality are encouraged to 
contact the corresponding author. The source code can be made 
available on request. 
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Figure 1 : Within Offline, tlie detector and event data structures are clearly sepai'ated. Both data structures have been complemented with the analogous classes for 
radio detection (marked in blue). In both the RDetector and the REvent, classes for Stations and Channels are present. Those in the detector data structures provide 
access to the detector description, the ones in the event data structures store data applying to specific radio events. 
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Figure 2: At both the Station and Channel levels of the REvent, data structures exist to store time-series and frequency-domain data. These are encapsulated in 
FFTDataContainers which transparently and efficiently handle all necessary FFTs without explicit interaction from the end-user. 
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Figure 7: Time traces (left) and frequency spectra (riglit) of tlie simulated event for the east and north channels. Using the simulated antenna characteristics 
(including the ampHfication by the LNA), the three-dimensional electric field vector has been projected to the two measurement channels. After the projection, 
white noise has been added by the RdChannelNoiseGenerator. The signal depicted here is what would be measured at the antenna foot-points over the whole 
frequency bandwidth. Please note that the frequency spectra correspond to the complete trace and not the zoomed-in time window shown here. 
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Figure 8: Time traces (left) and frequency spectra (right) of the simulated event for the east and north channels after folding in the channel response (amplifiers, 
filters, cable). Note that the cable delays have shifted the time pulses to later times. Also, the spectral bandwidth has been limited to the design bandwidth of 
the experimental channels, leading to a broadening of the pulses. The signal depicted here is what would be measured at the channel ADCs. Please note that the 
frequency spectra correspond to the complete trace and not the zoomed-in time window shown here. 
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Figure 9: Time traces (left) and frequency spectra (right) of the simulated event after all processing steps on the Channel level. The smoothing in the frequency 
spectra is due to the RdChannelRFISuppressor. These data again con'espond to the signal measured at the foot-points of the antennas, but this time limited to the 
40-70 MHz band. They are the starting point for the reconstruction of the three-dimensional electric field vector. Please note that the frequency spectra correspond 
to the complete trace and not the zoomed-in time window shown here. 
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Figure 10: Time traces (left) and frequency spectra (right) for the reconstructed three-dimensional electric field vector of the simulated event after application of the 
RdChannelWindowSetter. The leakage visible in the frequency spectra is due to the cutting of the time traces to a 500 ns window around the detected pulses. The 
frequency spectra con'espond to the time traces shown here. 
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Figure 1 1 : Time traces (left) and frequency spectra (right) of the simulated event for the final, Hann-windowed three-dimensional electric field vector. The significant 
vertical component of the raw simulation has been reconstructed by the analysis chain. The final reconstructed airival direction \s9 = 60.7° and (jy = 295.1°. The 
input direction for the simulated event was = 58.4° and ij) = 291.0°. The frequency spectra correspond to the time traces shown here. 



13 



